Is there something of the MCT in orientationally disordered crystals ? 
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Molecular Dynamics simulations have been performed on the orientationally disordered crystal 
chloroadamantane: a model system where dynamics are almost completely controlled by rotations. 
A critical temperature T c ~ 225 K as predicted by the Mode Coupling Theory can be clearly 
determined both in the a and (3 dynamical regimes. This investigation also shows the existence of 
a second remarkable dynamical crossover at the temperature T x > T c consistent with a previous 
NMR and MD study Jl| . This allows us to confirm clearly the existence of a 'landscape-influenced' 
regime occurring in the temperature range [T c — T x ] as recently proposed J2], |!|. 



In recent years, a great deal of work has been done to 
obtain a fundamental understanding of the glassy state 
and the glass formation. Despite these studies, no well- 
accepted theory has emerged so far At present, the 
only theory which gives precise predictions of the super- 
cooled liquids dynamics is the Mode Coupling Theory 
(MCT) |]. The idealized version of this theory particu- 
larly predicts an ergodic to nonergodic transition at the 
critical temperature T c which is usually depicted as an 
ideal-glass transition. However, experimentally, no sharp 
transition is observed at T c and the real glass formation 
only occurs at the calorimetric temperature T g < T c . 
It is suspected that some not so well understood pro- 
cesses i.e thermally activated hopping restore ergodic- 
ity below T c ('landscape-dominated' regime). During the 
past ten years, a variety of experiments |(| (time-of-flight 
neutron, Raman, depolarized light scattering, dielectric 
spectroscopy, optical Kerr effect) and Molecular Dynam- 
ics (MD) computer simulations have been made to val- 
idate MCT. It has been particularly shown that MCT 
gives an outstanding good description of simple liquids 
(binary Lennard- Jones, hard sphere) while for molecu- 
lar liquids where both translational and orientational de- 
grees of freedom (TDOF, ODOF) coexist, the predic- 
tions of the MCT seem only in fair agreement with the 
experimental or the simulation data. This raises partic- 
ularly the question, which remains a matter of debate at 
present, of the relative importance of rotations and trans- 
lations in molecular liquids during the dynamical slowing 
down. 

One class of compounds offers attractive possibilities 
to focus mainly on the role of the ODOF and provides 
a valuable alternative to the structural molecular liquid 
glass-formers. Indeed, some molecular crystals show a 
high temperature plastic phase in which the average po- 
sition of the centers of mass are ordered on a lattice 
while the orientations are dynamically disordered. Some 
of them, called glassy crystals 0, such as cyanoadaman- 
tane II , ethanol || [T(| or cyclooctanol 1 1 , can be deeply 
supercooled and present many properties characteristic 
of the conventional molecular liquid glasses such as a 



step in the specific heat at the glass transition or a non- 
Arrhenius behaviour of the relaxation times. 

In the present Letter, we discuss precisely, by means of 
a case study of chloroadamantane, the validity of MCT 
in orientationally disordered crystals. Indeed, although 
many experimental facts have been collected for this class 
of compounds in the deep supercooled domain, we have 
to emphasize that it really exist very few informations 
of dynamics far above T g . This particularly raises the 
question if MCT could be valid for such systems where 
dynamics are mainly controlled by rotations. One have 
to mention a numerical investigation made by Renner et 
at in |T3] where a system of colliding hard needles dis- 
tributed on a lattice, i.e a simple model of plastic crys- 
tal, has shown a rotational dynamical decoupling. Stim- 
ulated by this study and by the various predictions of 
MCT, we have chosen to study the chloroadamantane 
C10H15CI (noted Cla in the following) molecular crys- 
tal. This compound belongs to the substituted adaman- 
tane family which presents excellent experimental can- 
didates ||. Very recently, by means of NMR experi- 
ment and MD computer simulation, we have shown that 
this system exhibits a dynamical crossover transition at 
T x ~ 350 K in the pico- nanosecond regime jjj. It was 
interpreted as an indication of a a change of the en- 
ergy landscape topography as it has been recently pro- 
posed . In addition, this crossover was suspected to be 
the analogue of the Goldstein crossing temperature be- 
tween quasi-free diffusion and activated regime predicted 
in liquids. Two-step (a — (3) relaxations was also ob- 
served to emerge but there were some indications that 
T x probably did not correspond to the T c of the MCT as 
developed for structural glass-formers. 

MD simulation of orientationally disordered molecular 
crystal Cla have been performed at different tempera- 
tures from T = 220 to 500 K. Owing to the very long 
MD runs, equilibrium of the system could not be ob- 
tained below T ~ 220 K. A simple model was used which 
is completely described in |l^, [uj, so only the essential 
details are given in [ ^5| . 

Fig. [l] shows the orientational self intermediate scat- 
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FIG. 1: Orientational self intermediate scattering function 
F£(q,t) as function of time t for wave vector q — 3.06 A" 1 . 
F°(q,t) is denned as Eili Ef=i (exp[iq.(ui, a (t) - Ui,a(0))]) 
where Ui <a — Ri. a — Rf m . Rf m is the center of mass of 
molecule i and Ri, a denotes the position of site a in each 
molecules (N a = 2 in the Dresent studv). 
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FIG. 2: Elastic incoherent structure factor (EISF) as 
function of temperature for arbitrary wave vector q — 
3.06, 2.45, 1.22, 0.62 and 0.30 A" 1 . For a rotational diffusion 
model (solid lines), the EISF for one scattering atom moving 
on a sphere of radius a can be written: An(a) ~ jo(q.a) where 
jo is the 0-order spherical Bessel function JL6|. 



tering function F£(q,t), for arbitrary wave vector q = 
3.06A -1 (q = (q, 0,0)), as defined by Lewis et al. 
in 1 17 where only the ODOF are taken into ac- 
count. Clearly, lowering the temperature, as already ob- 



served for the orientational correlation functions Ci (t) = 
W Hi=i,N ( p l in @h the rotator phase of 

Cla exhibits a two-step relaxation. At intermediate 
times, F*(q,t) shows a plateau-like region which reveals 
the existence of an orientational cage effect i. e the rota- 
tional analogue of the translational cage effect observed 
in liquids. This transient regime is followed by a slow 
decay to a non zero value called elastic incoherent struc- 



ture factor (EISF) and noted A (q) — limt_>oo F£ (<Z> *) 
Jl6| which gives information on the time-average orien- 
tational geometry of the molecular motions. A (q) is 
displayed in the figure ^| for different wave vectors and 
exhibits an unusual temperature dependence which con- 
firms, as already seen in JjJ, the existence of a change 
of the rotational dynamics. Indeed, at high tempera- 
ture, an isotropic diffusion rotation model gives a good 
description of dynamics. Lowering the temperature, dis- 
crepancies of this model are clearly seen in the figure ^ 
and a jump-like between preferred orientations type of 
motion are to be assumed. 

In the following, we try to describe our data in the 
MCT framework. The idealized version |j| of this theory 
predicts a two-step relaxation scenario (fast j3, slow a) of 
all the time dependent correlators <fi(q,t). MCT partic- 
ularly states the following points: At short time, <f>(q, t) 
decays to a plateau value, the so-called nonergodicity pa- 
rameter classically noted /£. The dynamical regime as- 
sociated to this plateau is called (3. It is centered around 
the rescaling time t a which is given by t a — to|cr| _1 / 2a 
where io is a characteristic microscopic time and a is pro- 
portional to \T — T c \. Above T c , the late (3 regime or the 
early a relaxation is described by the following power law 
(going beyond first order): 



4>(q,t)=fZ-hW.(t/T) b + hV>.(t/T) 



21, 



(1) 



where the first two terms correspond to the classical 
von Schweidler law, the last term is a second order cor- 
rection and t — io|c|~ 7 with 7 = l/2a + 1/26. At 
long time, MCT also predicts that the previous fast 
regime is followed by a slow relaxation, called a, with 
the characteristic time r. Furthermore, parameters a 
and b are temperature and q independent and related 
via T 2 (l - a) /r(l - 2a) = T 2 (l + b) /T(l + 2b) where T(x) 
is the gamma function. In order to verify the validity 
of those scaling laws, we first tried to fix the exponent 
b using the factorization theorem in the (3 regime where 
hopping processes are not supposed to be dominant. The 
ratio i? (i) = (0(t)-0(ti))/(0(ta)-^(*i)), where (t ls t 2 ) 
are two different arbitrary times in the [3 regime, has 
been calculated at different temperatures and for differ- 
ent correlators <j) {Fu(l>t) f° r several wave vectors and 
Ci(t) for I = 1,2 as defined in M). Results are displayed 
in the figure ^. Clearly, it exists a time domain where 
the different correlation functions collapse onto a master 
curve in the (3 regime as it is predicted by MCT. Assum- 
ing the correlator function <\> departs from the plateau 
with a von Schweidler law, it is possible to show that the 
master curve depends only on the exponent b. Using a 
fitting procedure performed at several different temper- 
atures for all correlators mentioned previously, we ob- 
tained the best results for b — 0.79 which corresponds 
to a — 0.36 and 7 = 2.02. Then, fixing these values, 
we performed individual fit of the different correlators 
in the late (3 regime using the equation |l| with free pa- 
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FIG. 3: Ratio R^t) = (<j>(t) - <t>(h)) / \<t>{t 2 ) - 0(<i)) where 
(ti,t2) are two different arbitrary times in the (3 regime at 
T = 260 K. The following correlation functions have been 
used: <j>(t) = F*(q, t) for wave vectors q = 3.06, 2.45, 1.22, 0.62 
and 0.30 A -1 and Mt) = Ci(t) for I = 1,2 angular correla- 
tion functions (see 111). Master curve (solid line) obtained by 
the von Schweidler law with the exponent b = 0.79 is also 
displayed. 
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FIG. 4: Total prefactor [h^] 1 ^'" 1 as function of temperature 
for all investigated correlators (see Fig. ^). Fitting procedure 
reveals a linear dependence (solid lines) of the total prefactor 
which extrapolates to zero at the critical temperature T c ~ 
225 ± 8 K for all correlators. 



rameters / g c , = h^{q).T- b and = h^{q).T- 2b . 
Only the temperatures below 330 K where the plateau 
region is clearly defined were considered. The total pref- 
actor /jW is displayed in the figure [|. is not shown 
in the present Letter but it exhibits the same tempera- 
ture dependence. A temperature dependence is clearly 
found which extrapolates to zero at the critical temper- 
ature T c ~ 225 ± 8 K for the different correlators. A 
deviation to the linear evolution is observed close to T c 
which could be attributed to the occurence of hopping 
processes. The nonergodicity parameter /° show a rela- 
tively smooth evolution as function of temperature but 
we have no clear evidence of a cusp at T c due to the lack 
of data at low temperature. 

In the a regime, we defined the relaxation time r q 
as the time it takes to the different correlators to de- 
cay from 1 to 1/e. It is assumed that any characteristic 
times belonging to the a regime show asymptotically the 
same temperature dependence r q ~ t. The t\ and T2 re- 
laxation times are defined similarly for the orientational 
correlation functions Ci = \_2- For a direct comparison, 
orientational self intermediate scattering functions have 
been rescaled with the EISF as [F* (t)-A (q)]/[l-A (q)] 
in order to decay with time from 1 to zero. According 
to MCT, the a relaxation times r g ~ 1 ^ 7 (T) should yield 
straight lines intersecting the abscissa at T = T c . Fig. || 
shows that for all investigated correlators this prediction 
holds well over a relatively large temperature range. Ex- 
trapolation of the temperature dependence gives a crit- 
ical temperature of T c ~ 225 ± 8 K consistent with the 
/3-regime analysis. About 40 K above T c , discrepancies 
with the MCT power law prediction are seen, as already 
shown in the fast regime, and certainly due to hopping 
processes. Furthermore, the time-temperature superpo- 



sition principle has been also checked and we found that 
all correlators fall onto a master curve in the a time do- 
main when plotted as function of the rescaled time t/r q 
(not display in the present Letter). However, at high 
temperatures, close to the dynamical crossover previously 
observed at T x in O , a clear deviation of the a relaxation 
from the MCT prediction is found. At T > T x ~ 350 K, 
we have calculated that our data can be very well re- 
produced by a simple Arrhenius law. The T x crossing 
temperature alters slightly for the different correlators. 
At present, we have no clear interpretation of this fea- 
ture. Nevertheless, as it has been nicely demonstrated by 
Sastry et al. § | for a binary Lennard-Jones MD simu- 
lation, we obtained here for plastic crystal a strong indi- 
cation of the existence of a 'landscape-influenced' regime 
occurring in the temperature range [T c — T x \. The dy- 
namics should be 'landscape-dominated' only below T c . 

In conclusion, for the first time, this study reveals that 
dynamics of one orientationally disorder crystal can be 
relatively well described by the idealized version of the 
MCT. It is noteworthy to remark that such an agreement 
was expected for those systems where dynamics are al- 
most completely controlled by rotations, since some very 
recent extensions of MCT to molecular systems || [H| 
have particularly shown that the asymptotic predictions 
of the idealized version of MCT still hold. For Cla plas- 
tic crystal, the same critical temperature T c ~ 225 K 
can be extracted in the a and regime which strongly 
supports the MCT description. Furthermore, this inves- 
tigation confirms the existence of a second remarkable 
dynamical crossover at the temperature T x > T c consis- 
tent with previous calorimetric [^0| and NMR and MD 
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a relaxation time r q 1 ' 1 as function of temperature 



FIG. 5: 

for all investigated correlators (see Fig. |3|). The exponent 
7 = 2.02 has been extracted from the MCT analysis in the 
j3 regime. Solid lines indicate MCT-fit using the linear law 
r~ lh ~ (T - T ). Deviation from the MCT prediction is 



found at T > T x ~ 350 K where r q 
simple Arrhenius law. 



can be fitted with a 



studies |Q. This allows us to determine precisely the 
temperature range [T x — T c ] of the influenced-landscape' 
regime in this rotator system as defined recently by Sas- 
try et al. [| for a model glass forming liquid. Our 
results obtained for Cla call for new investigations: (i) 
It will be of great interest to study the q dependence of 
the different parameters /£, hq or r q which has not been 
performed in the present Letter for Cla due to the limited 
number of investigated wave vectors, (ii) Experimental 
and numerical extensions of the present study to other 
well documented orientationally disordered crystals such 
as cyclooctanol or ethanol are needed to fully validate 
the MCT picture, (iii) Theoretical studies are highly 
desirable in order to clarify the precise nature of the dy- 
namical crossover occuring at T x which remains unclear 
at present. 
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